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The  paper  presents  an  elementary  reaction  based  solid  oxide  fuel  cell  (SOFC)  model  coupled  with  anodic 
elementary  heterogeneous  reactions  and  electrochemical  charge  transfer  reactions  for  CO/CO2  fuel  based 
on  an  anode  supported  button  cell.  The  model  is  calibrated  and  validated  using  experimental  data 
obtained  for  various  C0/C02  fuel  compositions  at  750,  800  and  850  °C.  The  comparison  shows  that  the 
modeling  results  agree  well  with  the  experimental  data.  The  effects  of  operating  conditions  on  the  cell 
performance  and  the  detailed  species  concentration  distribution  are  predicted.  Then,  the  carbon  deposi¬ 
tion  on  the  SOFC  anode  with  C0/C02  fuel  is  experimentally  measured  and  simulated  using  the  elementary 
reaction  model.  The  results  indicate  that  lower  temperature  and  lower  operation  voltage  are  helpful  to 
reduce  the  possibilities  of  carbon  deposition  on  Ni  particle  surfaces. 
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1.  Introduction 

Solid  oxide  fuel  cells  (SOFCs)  are  known  for  their  fuel  flexibility 
and  tolerance  to  carbon  monoxide,  CO,  besides  their  high  efficiency, 
low  environmental  impacts  and  system  simplicity  for  stationary  or 
mobile  applications.  Syngas,  a  mixture  of  carbon  monoxide,  CO, 
and  hydrogen,  H2,  derived  from  steam  reforming  of  coal,  biomass, 
methane,  or  other  hydrocarbons,  can  be  used  in  an  SOFC  directly  in 
generating  electricity. 

In  order  to  optimize  the  anode  performance,  it  is  important 
to  investigate  the  comprehensive  electrochemical  oxidation  pro¬ 
cesses  within  SOFC  anode  integrating  with  experimental  data  and 
kinetic  mechanisms.  Since  it  is  difficult  to  clarify  internal  parameter 
distribution  and  complex  physical  phenomenon  in  fuel  cell  only  by 
the  expensive,  time-consuming  and  labor-intensive  experiments,  a 
validated  mechanistic  model  will  be  an  effective  method  for  inter¬ 
nal  phenomenon  predictions,  performance  evaluation  and  optimal 
design  of  the  SOFCs. 

With  respect  to  performance,  it  is  generally  agreed  that  the 
electrochemical  reaction  rate  of  CO  is  lower  than  that  of  H2  [1-3]. 
When  pure  CO  is  supplied  to  a  porous  anode,  Jiang  and  Virkar  [4] 
reported  that  the  maximum  power  density  is  about  40%  of  that 
obtained  using  pure  H2.  The  wide  variety  of  experimental  con¬ 
ditions  and  materials  in  CO  electrochemical  oxidation  result  in 
different  reaction  mechanisms  and  different  limiting  steps  being 
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proposed  by  different  researchers  [5,6],  and  related  mechanistic 
model  has  been  developed  [7,8].  However,  it  should  be  noted  that 
compared  with  the  systematic  experimental  and  theoretical  stud¬ 
ies  of  H2  electrochemical  oxidation  [9-17],  relative  fewer  studies 
have  been  performed  for  the  CO  electrochemical  oxidation  mech¬ 
anism,  especially  for  the  elementary  reaction  level  models  of  SOFC 
with  C0/C02  fuel. 

In  this  paper,  an  elementary  reaction  based  SOFC  model  coupled 
with  anodic  elementary  heterogeneous  reactions  and  electro¬ 
chemical  elementary  charge  transfer  reactions  for  C0/C02  fuel  is 
developed  based  on  an  anode  support  button  cell.  The  developed 
model  is  calibrated  and  validated  using  experimental  data  obtained 
on  a  button  cell  for  various  C0/C02  fuel  compositions  at  750,  800 
and  850  °C.  Then,  the  effects  of  operating  conditions  on  the  cell  per¬ 
formance  and  the  detailed  species  concentration  distribution  are 
predicted.  Finally,  the  carbon  deposition  in  the  SOFC  anode  with 
C0/C02  fuel  is  experimentally  measured  and  simulated  using  the 
elementary  reaction  model. 

2.  Experiment 

2.1.  Anode-supported  button  cell  structure  and  fabrication 

An  anode-supported  SOFC  button  cell  made  by  SICCAS  (Shang¬ 
hai  Institute  of  Ceramics  Chinese  Academy  of  Sciences)  was 
employed  in  this  study.  It  consisted  of  a  Ni/YSZ  anode  support 
layer  (680  pirn),  a  Ni/ScSZ  anode  active  interlayer  (15  pirn),  a  ScSZ 
electrolyte  layer  (20  pirn),  and  a  lanthanum  strontium  manganate 
(LSM)/ScSZ  cathode  layer  (15  pum).  The  diameter  of  cathode  layer 
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Nomenclature 

a[  pre-exponential  factor  in  sticking  coefficient 

expression 

A  pre-exponential  factor  (cm,  mol,  s) 

b[  temperature  exponent  in  sticking  coefficient 

expression 

c  concentration  (mol  m-3 ) 

df  activation  energy  in  sticking  coefficient  expression 

(J  mol-1 ) 

D  diffusion  coefficient  (m2  s-1 ) 

E  activation  energy  (J  mol-1 ) 

ECa  parameter  in  Eq.  (23 )  ( 1 30,000 J  mol-1 ) 

F  Faraday  constant  (96,485  C  mol-1 ) 

io  exchange  current  density  (A  m-2 ) 

*o,an  anodic  length  specific  current  density  (A  m-1 ) 

/F  volumetric  Faradaic  current  (A  m-3 ) 

k  reaction  rate  constant  (m,  mol,  s) 

I<  number  of  species 

Mk,  Mj  molecular  weight  (kg  mol-1 ) 

n  temperature  exponent 

ne  number  of  electrons  participating  in  the  reaction 

nt  total  number  of  particles  per  unit  volume 

N  number  of  the  reactions 

p  pressure  (Pa) 

P  whole  range  connection  probabilities  of  same  kind 

particles 

Q.  source  term  of  charge  balance  equations  (A  m-3 ) 

r  average  pore  size  (m) 

rel  mean  radius  of  the  electronic  conductor  particle  (m) 
R  gas  constant  (8.314  J  mol-1 1<-1 ) 

Rk  source  term  of  mass  balance  equations  (kg  m-3  s-1 ) 
s  molar  production  rate  (mol  m-2  s-1 ) 

S°  initial  sticking  coefficient 

Seff  effective  reaction  area  per  unit  volume  (m2  m-3 ) 
SNi  Ni  active  surface  area  per  unit  volume  (m2  m-3 ) 
Stpb  TPB  active  area  per  unit  volume  (m2  m-3 ) 

t  time  (s) 

T  temperature  (I<) 

V  voltage  (V) 

Vk,  Vj  diffusion  volume 

W  molecular  weight  of  the  gas-phase  species  in  stick¬ 
ing  coefficient  reaction  (kg  mol-1 ) 
xk  molar  fraction  of  gas-phase  species  k 

Z  coordination  number 

Greek  letters 

a  transfer  coefficient 

P  tuning  parameter  (£2_1  nrr2) 

y  parameter  modeling  the  rate  constant  from  sticking 

coefficient 

r  surface  sites  density  (mol  m-2 ) 

e  porosity 

ski  parameter  modeling  the  species  coerage 

r]  overpotential  (V) 

0  contact  angle  between  the  electronic  and  ionic  con¬ 

ductors  particles  (rad) 

0k  surface  coverage  of  species  k 

l±ki  parameter  modeling  the  species  coverage 

v  stoichiometric  coefficient 

o  conductivity  (S  m-1 ) 

r  tortuosity 

X  species  symbol 


Subscripts 

ac  anode  chamber 

act  active  layer 

an  anode 

ca  cathode 

cc  cathode  chamber 

ec  electrochemical  reactions 

el  electronic  conductor  particles 

elec  electronic 

electrolyte  electrolyte 

g  gas-phase  species 

i  reactions  index 

io  ionic  conductor  particles 

ion  ionic 

k  species  index 

Kn  Knudsen 

mole  molecular 

OCV  open  circuit  voltage 

ref  reference 

s  surface  species 

sp  support  layer 

Superscripts 

0  parameter  at  equilibrium  conditions 

bulk  bulk  phase 

eff  effective 

TPB  three  phase  boundary 


was  1.3  cm  and  diameters  of  other  layers  were  all  2.6  cm.  Before 
testing,  silver  paste  was  reticulated  on  the  anode  and  cathode  sur¬ 
face  by  screen-printing  for  current  collection. 

2.2.  Test  setup 

The  experimental  setup  can  be  shown  in  Fig.  1(a).  The  button 
cell  was  located  at  the  end  of  two  coaxial  alumina  tubes  and  was 
impacted  by  an  alumina  plate  which  was  strained  by  springs.  The 
Pt  mesh  was  used  as  the  cathode  current  collector  and  was  fixed  to 
the  porous  cathode  with  silver  paste  screen-printed  on  the  surface. 
The  oxidant  flowed  into  the  inner  tube  to  the  cathode  and  passed 
through  the  porous  Pt  mesh.  A  Ni  felt  (thickness  2  mm)  was  fixed  to 
the  anode  support  layer  with  silver  paste  to  collect  anode  current. 
The  fuel  was  induced  to  the  anode  surrounding  by  an  alumina  tube. 
Due  to  the  porous  structure  of  Pt  mesh,  the  anode  gas  reached  the 
anode  easily.  For  both  anode  and  cathode,  Pt  wires  were  used  as 
voltage  and  current  probes.  A  glass  ring  was  used  as  sealant  to  sep¬ 
arate  the  anode  gas  and  cathode  gas.  All  the  devices  were  enclosed 
in  one  quartz  tube  and  heated  by  a  furnace  to  keep  the  operation 
temperature.  Pure  H2  was  sent  to  the  chamber  for  1  h  to  fully  reduce 
the  anode  at  the  flow  rate  of  50  seem.  Then,  CO/CO2  mixtures  with 
different  composition  proportions  were  used  as  fuel. 

2.3.  Cell  performance  tests 

The  button  cell  performances  were  measured  experimentally 
at  three  temperatures  (750/800/850  °C)  with  different  ratios  of  CO 
and  C02,  as  listed  in  Table  1.  Oxygen  was  used  as  oxidant.  The  flow 
rates  of  fuel  and  oxidant  were  both  kept  at  300  seem. 

The  polarization  curves  were  measured  using  four-probe 
method  with  an  electrochemical  workstation  (IM6ex,  Zahner- 
Elektrik  GmbH,  Germany).  Electrochemical  impedance  spec¬ 
troscopy  (EIS)  was  performed  using  amplitude  of  10  mV  over  the 
frequency  from  0.1  Hz  to  100  kHz.  The  ohmic  resistance  of  but- 
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Fig.  1.  Illustration  of  the  (a)  experimental  setup  and  (b)  ID  model  geometry  of  SOFC  button  cell. 


ton  cell  was  estimated  from  the  high  frequency  intercept  of  the 
impedance  curve.  The  measurements  were  initiated  when  the  tem¬ 
perature  changes  30  min  later  until  the  system  was  stable.  After 
the  cell  tests,  the  microstructure  of  anode  was  characterized  using 
mercury  porosimeter. 


2.4.  Experimental  characterization  for  surface  carbon  deposition 

In  order  to  study  the  effects  of  temperature  on  the  anodic 
carbon  depositions,  a  mixture  of  C0/C02  with  90%  CO  and  10% 
C02  was  used  as  fuel  and  the  cell  was  kept  discharging  at  0.7  V, 
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Table  1 

Gas  composition  of  the  fuel  in  the  experimental  measurements. 


Temperature  (°C) 

Fuel  type 

Fuel  composition 

CO  (%) 

C02 (%) 

F750-1 

80 

10 

750 

F750-2 

40 

20 

F750-3 

20 

40 

F750-4 

10 

80 

F800-1 

80 

10 

800 

F800-2 

40 

20 

F800-3 

20 

40 

F800-4 

10 

80 

F850-1 

80 

10 

850 

F850-2 

40 

20 

F850-3 

20 

40 

F850-4 

10 

80 

750/800/850  °C,  for  2  h  after  the  cell  anode  is  fully  reduced  as  shown 
in  Table  2.  Then,  in  order  to  study  the  effects  of  fuel  compositions 
on  the  anodic  carbon  depositions,  fuel  with  different  C0/C02  ratios 
was  used  as  fuel  and  the  cell  was  kept  discharging  at  0.7  V,  750  °C, 
for  8h  after  the  cell  anode  is  fully  reduced.  It  should  be  noted 
that,  since  the  long  time  operation  at  higher  temperature  will  have 
effects  on  cell  structures  which  deteriorate  the  cell  performance. 
Thus,  we  choose  2  h  discharging  time  for  the  experiments  set  I  and 
8  h  discharging  time  for  the  experiments  set  II. 

After  the  carbon  deposition  operation,  the  cell  was  cooled  down 
to  ambient  with  Ar  (150  ml  min-1)  and  H2  (3  ml  min-1)  used  as 
protecting  gas.  About  2%  H2  was  added  to  maintain  a  reducing 
atmosphere  in  the  anode  chamber  to  prevent  oxidization  of  the 
deposited  carbon  within  the  anode. 

The  surfaces  of  the  anode  cross-sections  were  analyzed  by  a  X- 
ray  photoelectron  spectroscopy  (XPS)  (PHI  Quantera,  ULVAC-PHI, 
Kanagawa,  Japan)  with  a  monochromatic  Al  Ka  (hv  =  1486.7  eV) 
X-ray  source.  To  avoid  any  influence  of  carbon  contamination  dur¬ 
ing  the  sample  transfer  and  handing,  all  the  sample  surfaces  were 
peeled  1.2  nm  by  an  electron  gun  before  characterization.  The  spot 
size  for  analysis  was  300  [xm  x  300  [xm.  To  compensate  for  the 
surface-charge  effects,  the  binding  energy  scale  was  calibrated  with 
reference  to  the  binding  energy  of  0 1  s  at  530.2  eV.  The  survey  scans 
were  acquired  between  1200  and  0  eV.  The  elemental  surface  com¬ 
position  was  calculated  using  the  transmission  values  with  relative 
sensitivity  factors  specific  for  the  instrument  equipped  with  an  Al 
source. 

3.  Model  development 

3.1.  Model  assumptions  and  geometry 

The  following  model  assumptions  were  shown  as  follows: 


Table  2 

Control  experiments  for  the  anode  carbon  deposition. 

Experiments  Operating 

conditions 

Temperature  Fuel  com-  Discharging  Discharging 

positon  voltage  time 

Set  I  750/800/850  °C  90%  CO,  0.7  V  2h 

10%  C02 
90%  CO, 

Spf 

J  750°C  10%CO2  0.7  V  8  h 

50%  CO, 

50%  C02 
10%  CO, 

90%  C02 


(1)  Gases  are  assumed  as  ideal  gases. 

(2)  The  temperature  of  the  cell  is  uniform.  All  parameters  are  eval¬ 
uated  at  the  given  temperature. 

(3)  The  electrochemical  mechanism  is  modeled  using  a  set  of  ele¬ 
mentary  reactions  that  represent  chemical  reactivity  at  the 
molecular  scale.  Heterogeneous  thermochemical  and  electro¬ 
chemical  reactions  are  assumed  to  take  place  on  both  of  the  Ni 
surface  and  YSZ  surface.  One  adsorbed  CO  molecule  occupies 
one  vacant  position  on  the  surface,  neglecting  the  possibility  of 
two  vacant  positions  being 

(4)  The  charge  transfer  reaction  is  assumed  to  be  a  oxy¬ 
gen  spillover  reaction  taking  place  at  the  TPB  as: 
CO(Ni)  +  02-(YSZ)  <=►  C02(Ni)  +  (YSZ)  +  2e- 

(5)  The  microstructures  of  electrodes  are  stable  and  homogeneous 
in  operation.  The  carbon  deposition  influence  to  the  pore  struc¬ 
ture  and  reaction  activities  in  anode  is  disregarded.  Besides, 
the  distributions  of  the  two  conducting  phases  (electronic  and 
ionic)  in  electrodes  are  assumed  to  be  uniform. 

(6)  For  anode  heterogeneous  reactions,  mean  field  approximation 
is  used,  i.e.  the  surface  adsorbates  are  assumed  to  be  uni¬ 
formly  distributed  over  the  catalyst  surface.  And  the  transport 
of  surface  species  over  microscopic  distances  caused  by  surface 
diffusion  is  assumed  to  be  negligible. 

One-dimensional  (ID)  geometry  is  used  as  illustrated  in 
Fig.  1  (b).  The  model  calculation  domains  and  boundaries  are  labeled 
schematically  in  the  figure.  With  the  above  assumptions  and 
simplified  model  geometry,  a  ID  SOFC  model  is  developed  by 
considering  anodic  heterogeneous  chemistry,  charge  transfer  pro¬ 
cesses,  charge  balance  as  well  as  mass  balance. 

3.2.  Governing  equations 

3.2.1.  Anode  heterogeneous  chemistry 

Heterogeneous  chemistry  at  the  catalytic  surface  of  the  anode  is 
used.  Ni  is  an  effective  catalyst  for  surface  reactions,  especially  for 
hydrocarbon  fueled  SOFC.  A  heterogeneous  mechanism  simplified 
from  the  works  of  Hecht  et  al.  [18],  Janardhanan  and  Deutschmann 
[19]  and  Zhu  and  Kee  [20]  is  used,  as  shown  in  Table  3. 

The  surface  adsorbates  are  assumed  to  be  uniformly  distributed 
over  the  Ni  surface.  The  species  molar  production  rates  depend  on 
the  gaseous  species  concentrations  and  the  surface  species  concen¬ 
trations,  which  are  expressed  by  the  coverage.  The  coverage  0k  is 
the  fraction  of  the  surface  sites  covered  by  the  adsorbed  species  k.  It 
is  assumed  that  the  total  number  of  surface  active  sites  is  conserved 
and  the  saturation  sorbent  capacity  is  described  by  the  maximum 
surface  sites  density  r  [21  ].  The  uncovered  Ni  surface  is  treated  as 
a  dummy  surface  species. 

The  gaseous  adsorption-desorption  reactions  and  surface  reac¬ 
tions  are  written  in  the  general  form: 

I<g+Ks  Kg+I<s 

V'kXk  =>  VkXk  (1) 

k= 1  k=  1 

where  Xk  is  the  kth  species,  v'k  and  v"  are  the  stoichiometric  coef¬ 
ficients  of  the  reactants  and  products,  I<g  and  Ks  are  the  number 
of  gaseous  species  and  surface  species,  respectively.  The  net  molar 
production  rate  sk  of  a  gaseous  species  or  a  surface  species  in  a 
heterogeneous  reaction  are  written  as: 

N  Kg+Ks 

Sk  =  Y.(  vki  -  v’kOki  n  ckU  (2) 

1=1  k= 1 

where  N  is  the  total  number  of  reactions  and  ck  is  the  concentration 
of  the  kth  species.  For  all  the  surface  reactions  and  desorption  reac¬ 
tions,  the  reaction  rate  constant  Jq  for  the  ith  reaction  is  presented 
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Table  3 

Heterogeneous  reactions  mechanism  on  the  Ni  surface. 


Reaction 

Aa  (cm,  mol,  s) 

na 

Ea  (kj  mol-1) 

Adsorption 

If 

02  +  Ni(s)  +  Ni(s)^0(s)  +  0(s) 

1.000  x  10-  02b 

0.0 

0.00 

2f 

C02  +  Ni(s)  C02(s) 

1.000  x  10- 05b 

0.0 

0.00 

3f 

CO  +  Ni(s)^CO(s) 

5.000  x  10- 01b 

0.0 

0.00 

Desorption 

lb 

0(s)  +  0(s)^  Ni(s)  +  Ni(s)  +  02 

4.283  x  10  +  23 

0.0 

474.95 

2b 

C02(s)  -*  C02  +  Ni(s) 

6.447  x  10  +  07 

0.0 

25.98 

3b 

CO(s)^CO  +  Ni(s) 

3.563  x  10  +  11 
eco(s) 

0.0 

111.27 

-50.00c 

Surface  reactions 

4f 

C(s)  +  0(s)^CO(s)  +  Ni(s) 

5.200  x  10  +  23 

0.0 

148.10 

4b 

CO(s)  +  Ni(s)^C(s)  +  0(s) 

1.354x  10  +  22 

0CO(s) 

-3.0 

116.12 

-50.00c 

5f 

C0(s)  +  0(s)^C02(s)  +  Ni(s) 

2.000  x  10+19 
eco(s) 

0.0 

123.60 

-50.00c 

5b 

C02(s)  +  Ni(s)  CO(s)  +  O(s) 

4.653  x  10  +  23 

-1.0 

89.32 

a  Arrhenius  parameters  for  the  rate  constants  are  written  in  the  form:  k=ATnexp(-ElRT). 
b  Sticking  coefficient. 
c  Coverage-dependent  activation  energy. 


in  the  Arrhenius  form: 

Kg  -Ks  /  \ 

K  =  AjTni  exp  II  0) exp  J  (3) 

k=  1  '  ' 

where  A,-,  n{  and  E;  are  the  pre-exponential  factor,  temperature 
exponent  and  activation  energy  listed  in  Table  3,  R  is  the  gas 
constant,  T  is  the  temperature, ski  describes  the  species  coverage- 
dependency  of  the  rate  constant.  For  most  reactions,  which  are 
independent  of  species  coverage,  ski  is  zero.  For  reactions  3b,  4b 
and  5f,  the  reaction  rate  constants  depend  on  the  CO(s)  coverage, 

0CO(s)- 

For  adsorption  reactions,  the  rate  constants  are  expressed  in 
terms  of  the  sticking  coefficient  form: 


_  s°  /  RT 
~  Ty  y  2  tzW 


(4) 


where  S?  is  the  initial  sticking  coefficient  and  W  is  the  molecular 
weight  of  the  gas-phase  species.  And 


I<s 

y  =  Y^V'ki  (5) 

k 

Sticking  coefficient  S?  is  temperature  dependent  and  expressed  as 

S?  =  ajTbi  exp(-|f)  (6) 

where  a*  and  bi  are  dimensionless  parameters  and  d{  has  the  com¬ 
patible  units  with  RT.  These  parameters  can  be  treated  as  Arrhenius 
number  which  also  listed  in  Table  3. 

For  the  YSZ  surface,  few  validated  reaction  mechanisms  are 
available  in  the  literature.  Vogler  et  al.  [17]  considered  molecular 
adsorption  and  desorption  of  water,  water  dissociation,  and  bulk- 
surface  exchange.  In  this  study,  we  neglect  CO  and  C02  adsorption 
and  desorption  on  the  YSZ  surface,  and  only  consider  the  bulk- 
surface  exchange  of  oxygen.  Transport  of  bulk  oxygen  species  takes 
place  by  a  vacancy  diffusion  mechanism.  The  reaction  formula  and 
reaction  data  is  given  by 


02“(YSZ)  +  (YSZ,  bulk)  (YSZ)  +  02“(YSZ,  bulk) 

YSZ,  lb 


(7) 


where  kysz.if  is  determined  from  mass  action  kinetics, 
1<i  =  ki  0  exp(  -  ( Ei/RT ))  with  ki0  =  1.6  and  22  cm2  mol-1  s-1 ,  and 


activation  energy  Et  =  90, 900 J  mol-1,  kysz.ib  is  determined  from 
the  equilibrium  constant. 

3.2.2.  Anode  electrochemistry 

Theoretical  formulation  and  computational  modeling  of  ele¬ 
mentary  charge-transfer  chemistry  have  already  been  suggested 
from  different  researchers  [22,16,23-25].  Here,  the  anode  electro¬ 
chemical  reactions  were  taken  from  the  model  on  the  patterned 
SOFC  anode  [26]  and  validated  by  the  polarization  curves  and  EIS 
experimental  results.  CO  charge  transfer  pathway  is  considered  in 
this  study: 

„  kct,CO,f 

CO(Ni)  +  02_(YSZ)  ==  C02(Ni)  +  (YSZ)  +  2e“  (8) 

kct,CO,b 

where  /<ct,co,f  and  /<ct,co,b  are  the  forward  and  backward  reaction 
rates  CO  charge  transfer  reaction. 

According  to  Faraday’s  law,  the  current  source  at  the  anode,  and 
depending  on  the  charge  transfer  step,  can  be  expressed  as  follows: 

Qan  =  io,an^TPB  (9) 

k),an  =  2F(/<ct,f,COc02-(YSZ)CCO(Ni)  “  ^ct,b,COc(YSZ)cC02(Ni))  (1°) 

where  Qan  is  the  volumetric  specific  Faradic  current,  F  is  the  Fara¬ 
day  constant,  LTPB  is  TPB  length  per  unit  volume.  i0,an  is  the  length 
specific  current  density.  cC0(Ni)  and  Q:o2(Ni)  denote  the  surface  con¬ 
centrations  of  CO  and  C02  adsorbed  on  the  Ni  surface,  respectively. 
c02-(ysz)  and  c(ysz)  denote  volumetric  concentrations  of  oxygen 
interstitial  and  oxygen  vacancy  in  the  ionic  conductor. 

The  forward  and  backward  charge  transfer  reaction  rates  are 
written  as  follows: 


%u 

II 

u 

exp  ^ 

-(1  -oi)neF  \ 

RT  ^An) 

(11) 

:ct,b,i  -  fcct,b,i 

i  exp  | 

(cmeF  \ 

{  RT  11  an) 

(12) 

where  a  is  the  charge  transfer  coefficient.  Here,  it  should  be 
noted  that  the  global  representations  of  charge-transfer  chemistry 
the  symmetry  factors  do  not  sum  to  unit.  However,  to  simplify 
the  model  parameters,  the  single  charge  transfer  coefficient  was 
adopted  instead  of  anodic  and  cathodic  symmetry  factors.  r]a n  is 
the  anodic  overpotential,  which  is  defined  as: 

^7an  =  ^elec,an  —  ^ion,an  —  ^ref,an 


(13) 
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In  this  study,  the  anode  reference  potential  Vrefan  was  set  to 
zero.  The  Velecan  is  the  electronic  potential  of  nickel  at  the  TPB  inter¬ 
face  and  the  Vionan  is  the  ionic  potential  at  the  TPB  interface.  The 
parameters  kj?  f  co,  b  co  are  calculated  from: 


fc° 


*o,co 


ct,  f,  CO  " 


2H.TPBC°2_(Ysz)C°0(Ni) 


k°, 


*0,CO 


ct,b,CO“ 


2FLtpbc° 


(14) 


(YSZrC02(Ni) 


where  i0  is  the  exchange  current  density,  c°  denotes  the  species  sur¬ 
face  concentrations  at  equilibrium.  With  volumetric  concentrations 
of  oxygen  interstitial  and  oxygen  vacancy  in  the  ionic  conductor, 
the  equilibrium  surface  concentration  of  CO(Ni)  and  C02(Ni)  can 
be  determined  by  making  the  calculate  open  circuit  voltage  agree 
well  the  experimental  open  circuit  voltage  after  iterations. 

The  exchange  current  density  here  can  be  fitted  in  Arrhenius 
form: 


*’o,co  = 


AGcq\ 

RT  J 


(15) 


where  k  and  AG  are  the  tuning  parameters  in  Arrhenius  represen¬ 
tation  when  comparing  the  simulated  polarization  curve  with  the 
experimental  results. 


3.2.3.  Charge  balance 

For  the  electrodes,  the  charge  balance  equations  can  be  formu¬ 
lated  as 

v(— CTeffW,)=Q  (16) 

where  Q.  is  the  current  source,  V;  and  Vj  are  the  electric  potential  of 
either  ionic  or  electron  conductor  phases,  <reff  is  the  corresponding 
conductor  phase  effective  conductivity. 

The  governing  equations  for  the  electrode  charge  balances  are 
summarized  as  follows: 


•  Ionic  charge  at  the  cathode: 


v-(- 


wi01 


—  Qioi 


where  Ec a  is  1 30, 000 J  mol-1  and  /3  =  6.17  x  1011  m-2,  which 

indicate  the  oxygen  interface  conductivity  [27].  pc a  is  the  cathode 
local  overpotential,  which  is  defined  as: 

*7ca  =  ^elec.ca  —  ^ion.ca  —  Wef.ca  (22) 

where  Vrefca  is  the  cathode  local  relative  potential  difference 
between  the  electronic  conductors  at  a  reference  state.  Similar  to 
the  cathode  local  overpotential,  the  anode  local  overpotential  pan 
is  defined  as 

?7an  =  ^elec.an  —  ^ion.an  —  ^ref.an  (23) 

By  setting  anode  reference  potential  Vrefan  to  zero,  the  cathode 
reference  potential  Vref  ca  equals  to  the  actual  cell  open  circuit  volt¬ 
age  (OCV)  Vocv-  Essentially,  the  SOFC  is  a  concentration  cell  of  02 
and  the  OCV  can  be  determined  by  Nernst  equation  as  [28]: 


where  p§ }  and  p°2  are  the  equilibrium  oxygen  partial  pressures 
in  cathode  and  anode,  respectively.  pa  2  can  be  determined  by  the 
heterogeneous  reactions  and  the  compositions  of  fuel  in  anode. 

The  electrolyte  charge  balance  equation  can  be  briefly  expressed 
as 

V  '  (_afon,electrolyteV^ion, electrolyte)  =  0  (25) 

where  cr^fn  electrolyte  is  the  effective  ionic  conductivity  of  electrolyte 
and  Vion  electrolyte  is  the  ionic  potential  in  electrolyte. 

3.2.4.  Mass  balance 

The  mass  balance  equation  in  porous  electrode  can  be  formu¬ 
lated  as  [29]: 

V(-DfVck,g)  =%g  (26) 

where  s  is  the  porosity  of  electrode,  ckg  is  the  gas  molar  concentra¬ 
tion,  Rkg  is  the  mass  balance  source  term  of  gaseous  specie  within 
the  porous  electrode,  and  is  the  effective  diffusivity  of  gaseous 
species  /<,  which  can  be  formulated  as 


•  Electronic  charge  at  the  cathode: 

^  '  (— ^elec.ca^^^c.c3^  =  Qelec.ca  =  —  Qion.ca  (18) 

•  Ionic  charge  at  the  anode: 

V  •  (-CTfon,anVVVan)  =  Qion.an  =  Qan  (19) 

•  Electronic  charge  at  the  anode: 

^  '  (— ^elec.an^^elec.an)  —  Qelec.an  —  —  Qan  (20) 

The  Qan  is  the  volumetric  Faradic  current  which  can  be  found 
in  Section  3.2.2.  c™  and  cj™lk  are  the  cathode  oxygen  concen¬ 
trations  at  the  TPB  and  in  the  bulk,  respectively,  a  is  the  charge 
transfer  coefficient,  ne  is  the  number  of  electrons  participating  in 
the  reaction. 

The  cathode  exchange  current  density  io,Ca  is  expressed  as 
io,ca=^  exp(-^)(p§)0'25  (21) 


Df  = 


'  +  ' 


Deff  Deff 

,uk,  mole  uk,Kn 


(27) 


where  }  and  £>jf[<n  are  the  effective  molecular  diffusion  coef¬ 
ficient  and  effective  Knudsen  diffusion  coefficient,  respectively.  For 
multi-components  gases,  the  effective  molecular  diffusion  coeffi¬ 
cient  can  be  calculated  as: 


D 


eff 

/c,mole 


O -**)/>  ,  ^ 


j  =  i 


K 


(28) 


where  xk  and  Xj  are  the  fractions  of  molar  concentrations  of  gaseous 
species  k  and  j,  and  D ^  is  the  effective  binary  molecular  diffu¬ 
sion  coefficient.  The  effective  molecular  diffusion  coefficient  of  kth 
gaseous  species  is  restricted  to  the  situation  wherein  the  species  k 
diffuses  in  a  mixture  of  stagnant  and  non-transferring  species.  Con¬ 
sidering  the  porous  material  property,  effective  binary  molecular 
diffusion  coefficient  Dfff  and  effective  Knudsen  diffusion  coefficient 

D£f[<n  can  be  calculated  as 


Df:=b~DkJ  = 


k,j 


0.00101eT1-75((l/Mfc)  +  (1/Mj))1/2 


zp 


vf3  +  v)/3 


(29) 
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Table  4 

Boundary  conditions. 


Boundary 

Ionic  charge 

Electronic 

charge 

Mass 

balance 

d^cc/fuel 

Insulation 

0 

Cg,an  for 
gaseous 
species 

9^2an.sp/cc 

Insulation 

Continuity 

Insulation 
for  surface 
species 

^^ari-act/an-sp 

Continuity 

Continuity 

Continuity 

9  -^electro  lyte/an.act 

Continuity 

Insulation 

Insulation 

9i£2Ca/eIectrolyte 

Continuity 

Insulation 

Insulation 

9^2ca/air 

Insulation 

Vca 

Cg,ca 

neff  8  n 
Uk,  Kn  =  ~Uk,Kn  = 


4 6.  I  8 RT 
3  xV  y  7zMk 


(30) 


where  r  is  the  tortuosity  factor,  V  is  the  diffusion  volume.  The  value 
of  diffusion  volume  18.0  m3  mol-1  for  CO  and  26.7  m3  mol-1  for 
CO2.  M  is  the  molecular  weight,  p  is  the  total  pressure  of  gases,  and 
r  denotes  the  average  pore  radius. 

In  the  cathode,  the  relationship  between  the  mass  balance 
source  term  and  the  current  source  term  can  be  determined  by 
Faraday’s  law  as 


ro2  = 


Qelec.ca 

4F 


(31) 


The  reaction  rate  of  nitrogen  is  zero  and  the  molar  fraction  can 
be  determined  by 


xN2=1-Xo2  (32) 

Different  from  the  cathode,  the  mass  balance  source  terms  in 
the  anode  are  determined  by  kinetic  reaction  rates  per  unit  volume 
of  heterogeneous  and  electrochemical  reactions.  For  both  gaseous 
and  surface  species  in  the  anode,  the  source  terms  Rkg  and  Rk  s  can 
be  expressed  in  the  same  form: 

N  Kg+Ks 

Rk  =  S'*  ■  sk  =  S*Y,  K-  -  v'ki)  Ki  n  ckU  (33) 

1=1  k=  1 

where  Seff  is  the  effective  reaction  area  per  unit  volume,  which 
equals  to  SNi  for  heterogeneous  reactions. 

In  addition,  since  the  concentrations  of  oxygen  interstitial  and 
oxygen  vacancy  in  ionic  conductor  are  several  orders  of  magnitude 
larger  than  the  concentrations  of  gases  and  surface  species,  the  vari¬ 
ation  of  concentrations  of  oxygen  interstitial  and  oxygen  vacancy 
is  neglected. 


3.3.  Boundary  conditions 


TPB  length  are  needed  in  the  calculation.  The  parameter  5Tpb  can 
be  formulated  by  using  the  particle  coordination  number  in  binary 
random  packing  of  spheres  together  with  percolation  theory  as 
[30,31] 

Stpb  =  n  sin2  6»r2  NtotndnioZclZioPclPio/Z  (34) 

where  Z  is  the  mean  coordination  number,  rel  is  the  mean  radius 
of  the  electronic  conductor  particle,  0  is  the  contact  angle  between 
the  electronic  and  ionic  conductors  particles,  Ntot  is  the  total  num¬ 
ber  of  particles  per  unit  volume,  ne\  and  nl0  are  the  fraction  number 
of  electronic  and  ionic  conductor  particles,  Zel  and  Zio  are  the  coor¬ 
dination  numbers  of  electron  and  ion  conductor  particles,  and  Pel 
and  Pio  are  the  whole  range  connection  probabilities  of  the  same 
kind  particles.  Flere  Stpb  is  treated  as  the  contact  area  of  different 
types  of  particles.  The  detailed  model  description  and  parameter 
calculation  can  be  found  in  the  literature  [30,32]. 

Similar  to  the  calculation  Stpb,  the  parameter  LTpb  (TPB  length 
per  unit  volume)  can  be  calculated  as 

Ltpb  =  27rrel  sin 0Ntotne\Zel_ioPelPio  (35) 

The  parameter  SNi  (effective  Ni  surface  area  per  unit  volume) 
can  be  calculated  by  the  total  Ni  particles  surface  area  subtracting 
the  contact  surface  area  between  Ni  and  YSZ  particles  and  between 
two  Ni  particles  around  the  Ni  particles  per  unit  volume.  Then,  the 
effective  Ni  surface  area  per  unit  volume  can  be  formulated  as 

SNi  =  Trr^ntHe!  (4  -  sin2  PnioZelZio/Z  -  sin2  0nelZelZel/Z)  (36) 

The  pore  structure  of  the  anode  support  layer  was  charac¬ 
terized  using  mercury  porosimeter.  The  mean  pore  diameter, 
porosity  and  total  pore  area  were  found  to  be  0.387  [xm,  0.335 
and  8.54  x  106  m2  m-3,  respectively.  To  simplify  the  calculation, 
the  mean  particle  diameters  of  the  two  conductors  are  assumed 
to  be  the  same  and  equal  to  the  mean  pore  diameter  [32].  With 
this  assumption  and  the  expressions  of  effective  reaction  areas,  the 
calculated  values  of  Stpb  and  SNi  in  the  anode  support  layer  are 
2.22  x  105  and  3.97  x  106  m2  m-3,  respectively. 

It  is  very  difficult  to  characterize  the  pore  structures  of  anode 
active  layer  and  cathode  layer  by  experiment  since  they  are  very 
thin  and  hard  to  be  separated  from  the  button  cell.  We  have  used 
image  processing  software  to  determine  the  pore  size  and  porosity 
of  each  layer  compared  to  anode  support  layer  from  SEM  image 
based  on  quantitative  stereology.  It  is  found  that  the  mean  pore 
size  in  anode  support  layer  is  nearly  1.2  and  1.5  times  of  that  in 
cathode  layer  and  anode  active  layer,  but  the  porosities  of  all  three 
layers  are  almost  the  same.  Table  5  lists  the  porosity,  mean  pore 
diameter,  5Tpb  and  SNi  in  each  layer. 

Then,  the  material  conductivities  and  some  other  calculation 
parameters  can  be  found  in  Table  6. 


According  to  the  operation  conditions  and  model  simplifica¬ 
tions,  the  boundary  conditions  of  charge  and  mass  balances  partial 
differential  equations  are  listed  as  in  Table  4. 

The  boundary  conditions  “insulation”  and  “continuity”  mean 
that  the  partial  derivative  is  zero  and  the  flux  is  continuous  of 
the  variables  at  the  boundary.  cgian  and  cgiCa  in  the  table  are  the 
molar  fractions  of  gaseous  species  in  the  anode  and  cathode.  At 
the  interface  of  anode  and  anode  chamber,  the  boundary  condi¬ 
tions  for  surface  species  are  “insulation”,  which  is  different  from 
gaseous  species.  That  is  because  there  is  no  surface  species  out  of 
the  anode.  The  VCa  is  the  cell  operation  voltage  in  the  calculation. 


3.5.  Solution  method 

The  calculations  of  model  were  performed  using  the  finite  ele¬ 
ment  commercial  software  COMSOL  MULTIPHYSICS®,  Version  3.2. 

The  button  cell  performance  was  calculated  at  a  given  cell  volt¬ 
age  Vca.  For  ID  SOFC  model,  the  average  current  density  at  a  given 
cell  voltage  was  achieved  as  that  in  the  electrolyte  layer.  By  set¬ 
ting  different  cell  voltages,  a  complete  polarization  curve  can  be 
generated. 

4.  Results  and  discussion 


3.4.  Model  parameters 

From  the  model  development  section  of  charge  balance  and 
mass  balance  equations,  the  effective  reaction  area,  TPB  area  and 


4A.  Model  calibration  and  validation 

Since  some  model  parameters  cannot  be  determined  directly 
from  experiments  or  published  literature,  a  calibration  process  in 
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Table  5 

Pore  structure  parameters  in  porous  electrode. 


Cell  layer  Porosity  Mean  pore  diameter  (|jim) 

Stpb  (m2  m“3) 

Ltpb  (mm  3 ) 

SNi  (m2  m“3) 

Anode  support  layer  0.335 

0.193 

2.22  x  105 

2.76  x  107 

3.97  x  106 

Anode  active  layer  0.335 

0.129 

3.33  xlO5 

6.20  x  107 

5.96  xlO6 

Cathode  layer  0.335 

0.161 

2.66  xlO5 

3.97  x  107 

- 

Table  6 

Properties  and  parameters  for  model  calculation. 

Property  and  parameter 

Value  or  expression 

Unit 

Ionic  conductivity  (rrion) 

ScSZ 

6.92E4  exp(-9681  /7) 

Sm-1 

YSZ 

Electronic  conductivity  (crelec) 

3.34E4  exp(  - 1 0300/7) 

Sm-1 

LSM 

4.2E7  IT  exp(  - 1 1 50/7) 

Sm-1 

Ni 

3.27E6— 1065.37 

Sm-1 

Equivalent  ionic  conductivity  of  electrolyte  layer  (creiectroiyte) 

-3.622E-572  +  0.0837-  46.343a 

Sm-1 

Concentration  of  oxygen  interstitial  in  the  YSZ  (c02-(YSZ)) 

4.45E4 

mol  m2 

Concentration  of  oxygen  interstitial  in  the  YSZ  (c(Ysz)) 

4.65E3 

mol  m2 

Maximum  surface  sites  density  (T”) 

2.6E-5 

mol  m2 

Cathode  tortuosity  (rca) 

3.0 

Ni  felt  current  collector  porosity 

0.6b 

a  Experimentally  measured. 
b  As  received. 

a  reasonable  range  to  fit  the  experimental  data  is  necessary.  In  this 

the  OCV  and  cell  performance  decreases  significantly.  It  is  clear 

paper,  the  model  parameters  were  tuned  according  to  the  exper¬ 
imental  polarization  curves  with  fuel  F800-2  (40%  CO,  60%  C02) 
and  F800-4  (10%  CO,  90%  C02)  cases  at  800  °C.  The  tuning  param¬ 
eter  values  are  listed  in  Table  7.  Once  the  model  parameters  are 
determined,  they  were  not  changed  in  other  calculations. 

It  can  be  seen  that  the  tortuosity  of  both  the  anode  and  cathode  is 
set  at  3  and  the  tortuosity  of  the  Nickel  felt  layer  is  set  at  2.  It  seems 
reasonable  since  the  observed  range  for  porous  sintered  ceramics  is 
usually  2-10,  and  most  often  in  the  range  of  2-6  [33].  The  transfer 
coefficient  a  depends  on  the  symmetry  of  the  activation  barrier, 
usually  ranges  from  about  0.3  to  0.6. 

Fig.  2  shows  the  modeling  and  experimental  polarization  curves 
for  various  CO/C02  fuel  mixtures  at  750/800/850  °C.  It  can  be  seen 
that  the  modeling  polarization  curves  agree  well  with  the  exper¬ 
imental  data  at  800  °C  with  composition  10%  CO,  90%  C02  and 
composition  40%  CO,  60%  C02,  which  are  the  base  cases  for  model 
calibration.  It  should  be  noted  that  there  are  the  distinct  kinks 
around  0.5  V  in  each  polarization  curve.  However,  this  phenomenon 
is  not  considered  in  the  proposed  model.  As  a  preliminary  inter¬ 
pretation,  this  phenomenon  may  be  due  to  the  unstable  periodic 
reaction  rate  at  certain  operating  voltages  which  has  been  sug¬ 
gested  by  Holtappels  et  al.  [34].  As  the  CO  concentration  decreases, 


Table  7 

Model  tuning  parameters. 


Parameters 

Value 

Cathode  tortuosity 

3 

Cathode  electrochemical 
kinetics  parameter,  ft 
(ft-1  m-2) 

1.26  x  109 

Cathode  charge  transfer 
coefficient  (ofca) 

0.33 

Anode  tortuosity 

3 

Anode  current  collector  Ni  felt 
tortuosity 

2 

Anode  charge  transfer 
coefficient  (aran) 

0.6 

Anodic  electrochemical 
kinetics  parameter,  kco  (Am-1) 

1.68E15 

Anodic  electrochemical 
kinetics  parameter,  AGCo 
(Jmol-1) 

165615 

that  there  is  a  limiting  current  density  at  zero  voltage  caused  by 
anodic  concentration  polarization.  The  simulated  limiting  current 
densities  also  match  well  with  the  values  found  experimentally. 
While,  it  can  be  found  that  the  curvature  of  the  polarization  curves 
is  more  significant  with  the  decreasing  of  CO  content  due  to  the 
concentration  effects  at  lower  operating  voltage. 

The  polarization  curve  at  750  and  850  °C  has  similar  features 
to  that  at  800  °C.  At  750  and  850  °C,  the  modeling  results  agree 
well  with  the  experimental  data  at  low  current  density  but  devi¬ 
ate  from  experimental  data  for  some  cases  at  high  current  density, 
which  shown  in  Fig.  2(a)  and  (c).  The  modeling  results  at  750  °C 
underestimate  the  diffusion  resistance  of  concentration  polariza¬ 
tion  slightly  at  high  current  density.  For  higher  CO  content  (e.g. 
>50%)  at  750  °C,  the  concentration  polarization  is  not  significant. 
The  calculated  limiting  current  densities  are  lower  than  experimen¬ 
tal  data  at  850  °C,  which  is  more  pronounced  in  high  CO  content 
cases.  The  model  overestimates  the  diffusion  resistance  of  con¬ 
centration  polarization,  especially  at  higher  current  density.  The 
simplification  of  model  geometry  might  be  the  important  source  of 
the  derivation.  The  anode  of  button  cell  is  larger  than  the  cathode, 
which  result  in  the  non-uniform  distributions  of  gas  concentra¬ 
tion  in  anode  in  radial  direction.  In  fact,  the  model  is  not  exactly 
accorded  with  ID  geometry  assumption  and  this  assumption  will 
lead  to  the  deviation  in  calculation,  especially  for  the  high  current 
cases.  Comparison  between  the  modeling  results  and  experimental 
data  in  Fig.  2(a)-(c)  indicates  that  the  cell  performance  increases 
and  concentration  polarization  becomes  more  significant  with  the 
increasing  of  temperature  for  the  same  CO/C02  composition.  There¬ 
fore,  the  influence  of  geometry  simplification  is  more  obviously  at 
850  °C,  which  shown  in  Fig.  2(c).  In  addition,  the  neglect  of  bulk 
diffusion  in  anode  chamber  also  adds  contribution  to  the  devia¬ 
tion  between  the  calculated  results  and  experimental  data  at  high 
current  density. 


4.2.  Species  concentration  distributions  within  anode 

Since  the  electrode  is  very  thin  and  works  at  high  tempera¬ 
ture,  it  is  very  difficult  to  characterize  the  important  parameter 
distribution,  especially  hard  to  get  the  surface  reaction  elements 
information.  In  this  section,  the  elementary  surface  concentra- 
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Current  density(Arrf2) 


Fig.  2.  Modeling  and  experimental  polarization  curves  for  various  CO/CO2  fuel  mixtures  (a)  850  °C;  (b)  800  °C;  and  (c)  750  °C. 


tion  indicates  the  concentration  on  the  Ni  surface,  which  does  not 
include  the  distributed  information  of  elementary  species  near  the 
TPB. 

Fig.  3  shows  the  effects  of  temperature  on  the  surface  coverage 
of  NI(s),  CO(s),  H(s)  and  0(s)  for  fuel  mixture  80%  CO  and  20%  C02 
at  800  °C  and  cell  voltage  at  0.8  V.  It  can  be  seen  that  the  CO(Ni)  and 
(Ni)  (Ni  vacancies)  are  the  major  species  in  anode,  while  C02(Ni), 
O(Ni)  are  the  minor  species  which  is  at  the  order  of  le-5  to  le-6. 
The  surface  species  turn  to  equilibrium  quickly  near  the  anode  outer 
surface  as  gases  species  due  to  the  rapid  heterogeneous  reactions. 
The  calculation  results  indicated  that  the  coverage  of  CO(Ni)  and 
O(Ni)  decreases;  the  coverage  of  C02(Ni)  and  (Ni)  increases  with  the 
increasing  of  the  temperature  from  750  °C  to  850  °C.  In  addition,  it 
can  be  found  that  the  CO(Ni)  coverage  decreases  along  the  electrode 
thickness  direction.  In  fact,  the  electrochemical  reaction  goes  fast 
near  the  electrode/electrolyte  interface  and  thus  consumed  more 
CO(Ni)  surface  species  to  react  with  oxygen  ions. 

Fig.  4  shows  the  coverage  distribution  of  various  surface  species 
in  the  anode  for  fuel  mixture  80%  CO  and  20%  C02  at  800  °C  and 
cell  voltage  at  1.0,  0.8  and  0.6  V.  The  calculation  results  indicated 
that  variation  of  the  surface  species  distribution  along  the  electrode 
became  more  significant  at  lower  operating  voltage  due  to  the  rel¬ 
atively  faster  electrochemical  reaction  rates.  Especially  when  the 
cell  voltage  is  kept  at  0.6  V,  the  coverage  ratio  of  the  surface  species 
changes  rapidly.  This  indicated  that  the  species  surface  concentra¬ 
tion  is  much  larger  at  high  operating  mode  compared  with  that 


at  low  overpotential  operating  mode,  which  may  also  lead  to  the 
surface  diffusion  phenomenon  at  local  catalyst  surface. 

Fig.  5  further  shows  the  effects  of  fuel  compositions  on  the 
surface  coverage  of  (Ni),  CO(Ni),  O(Ni)  and  C02(Ni)  with  operat¬ 
ing  temperature  at  800  °C  and  cell  voltage  at  0.8  V.  It  can  be  seen 
that  the  surface  species  coverage  is  significantly  affected  by  the 
fuel  composition.  When  CO  content  increases  from  20%  to  80%,  the 
CO(Ni)  surface  coverage  get  larger,  while  the  (Ni),  C02(Ni)  and  O(Ni) 
contents  get  smaller. 

4.3.  Modeling  analysis  on  carbon  deposition  phenomenon 

The  effects  of  carbon  deposition  process  will  be  discussed  in 
both  experimental  and  modeling  aspects  in  this  section.  Flowever, 
it  should  be  noted  that  in  the  model,  only  the  carbon  coverage  on  the 
catalyst  surface  is  considered,  while  the  phenomenon  of  deposited 
carbon  growth  on  the  anode  surface  is  not  considered.  Thus,  the 
modeling  analysis  is  focused  on  the  qualitative  or  semi-qualitative 
but  not  quantitative.  In  addition,  since  it  is  hard  to  detect  the  car¬ 
bon  deposition  process  directly  from  experiments,  the  modeling 
analysis  will  be  significant  important  for  operating  strategy  choice 
and  for  understanding  the  detailed  electro-thermal  chemistry  with 
SOFC  anode. 

The  effect  of  temperature  on  carbon  deposition  degree  is  shown 
in  Table  8.  It  suggests  that  the  carbon  atomic  percentage  decreases 
with  the  temperature.  Then,  the  modeling  results  of  carbon  cov- 
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Distance  from  anode  surface  (  pm) 


Fig.  3.  Coverage  distributions  of  surface  species  within  anode  at  750/800/850  °C. 
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Fig.  4.  Coverage  distributions  of  surface  species  within  anode  at  different  operating 
voltages  (1.0V/0.8V/0.6V). 


Fig.  5.  Coverage  distributions  of  surface  species  within  anode  with  different  fuel 
compositions. 


Table  8 

Effects  of  temperature  on  anode  carbon  deposition  degree. 


Temperature  (°C) 

Atomic  percentage  (at.%)  within  anode 

C 

O 

Ni 

Zr 

Y 

750 

12.87 

56.56 

11.71 

16.54 

2.32 

800 

11.60 

58.18 

11.62 

16.38 

2.22 

850 

10.72 

57.52 

13.30 

16.39 

2.07 

erage  variation  at  different  temperatures  are  shown  in  Fig.  6. 
The  modeling  results  also  indicate  that  the  surface  coverage  of 
C(Ni)  increases  significantly  when  the  temperature  deceases  at 
certain  current  density,  which  is  a  good  qualitative  agreement 
with  the  experimental  results.  And,  it  also  suggests  that  the  car¬ 
bon  deposition  is  mitigated  at  higher  current  density.  In  fact,  at 
a  higher  current  density,  more  surface  specie  O(s)  is  generated 
in  electrochemical  reactions  and  hence  reducing  the  possibility 
of  carbon  deposition  on  Ni  surfaces.  Thus,  it  can  be  deduced  that 
the  carbon  deposition  degree  of  syngas  fueled  SOFC  could  be 
reduced  by  raising  the  operation  temperature  and  reducing  cell 
voltage. 

The  effects  of  fuel  composition  on  carbon  deposition  degree  at 
750  °C  are  shown  in  Table  9.  It  suggests  that  the  carbon  atomic  per- 


Table  9 

Effects  of  CO  mole  fraction  in  the  fuel  on  the  anode  carbon  deposition  degree. 


CO  mole  fraction  (%) 

Atomic  percentage  (at.%)  within  anode 

C 

O 

Ni 

Zr 

Y 

90 

14.32 

54.18 

12.78 

16.60 

2.13 

50 

12.82 

58.15 

10.28 

16.55 

2.21 

10 

10.77 

56.29 

12.18 

18.26 

2.51 
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Fig.  6.  Surface  carbon  species  converge  within  anode  at  750/800/850  °C. 
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Fig.  7.  Surface  carbon  species  converge  within  anode  at  different  CO  contents. 


centage  increases  with  the  mole  fraction  of  CO  in  the  fuel.  Fig.  7 
further  shows  the  modeling  results  of  fuel  composition  effects  on 
carbon  deposition  degree.  It  suggests  that  the  degree  of  carbon 
deposition  aggravates  as  the  CO  content  decreases.  This  qualitative 
trend  agreed  well  with  the  experimental  observations.  It  also  indi¬ 
cated  that  when  the  CO  mole  fraction  is  less  than  50%,  the  C(Ni) 
coverage  ratio  can  be  kept  at  the  level  under  1.5e-4.  However, 
when  the  CO  mole  fraction  increased  to  90%,  the  C(Ni)  coverage 
ratio  increases  fast  to  le-3  level,  which  is  almost  6  times  larger 
than  that  at  50%  CO  fuel  composition.  It  further  indicated  the  neces¬ 
sity  for  keeping  the  CO  content  at  a  lower  level  to  release  the  carbon 
deposition  risk. 

5.  Conclusion 

A  detailed  one-dimension  elementary  reaction  model  of  an 
anode-supported  SOFC  operating  with  CO/CO2  mixtures  based 
on  button  cell  geometry  was  developed.  The  model  incorpo¬ 
rated  anodic  elementary  heterogeneous  reactions,  electrochemical 
kinetics,  electrodes  microstructure  and  complex  transport  phe¬ 
nomena  (mass  and  charge  transport)  within  electrode. 

The  model  was  calibrated  and  validated  with  the  experimen¬ 
tal  data  with  different  CO  and  CO  compositions  at  750,  800  and 


850  °C.  The  results  indicated  that  the  model  agreed  reasonably  well 
with  the  experimental  data.  The  model  can  well  be  applied  for 
the  detailed  understanding  of  the  electro-thermal  chemical  and 
transport  processes  with  porous  electrode  and  on  electrode  par¬ 
ticles  surface,  and  also  the  model  can  be  used  in  the  cell  operating 
condition  optimization. 

For  the  heterogeneous  reactions,  CO(Ni)  and  Ni  vacancies  are 
the  major  species  in  anode,  while  C(Ni),  C02(Ni)  and  O(Ni)  are 
the  minor  species  with  the  coverage  is  about  10-5  to  10-6.  As 
the  cell  voltage  decreases,  the  faster  electrochemical  reactions 
resulting  in  the  increasing  of  C02(Ni)and  O(Ni)  and  the  decreas¬ 
ing  of  CO(Ni),  C(Ni).  The  effects  became  more  obvious  near  the 
TPB.  The  coverage  of  CO(s)  and  O(s)  decreases  with  increasing 
temperature  due  to  the  higher  desorption  rates.  In  addition,  both 
of  the  modeling  and  experimental  results  indicate  that  lower 
temperature  and  operation  voltage  are  helpful  to  reduce  the  pos¬ 
sibility  of  carbon  deposition  on  Ni  surfaces.  While,  as  the  C02 
content  in  the  fuel  increases,  the  degree  of  carbon  deposition 
aggravates. 
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